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Abstract 



The normalized state \ip(t)) = ci(i)|l) + C2 (i) |2) of a single two-level system performs 
\J2 ', oscillations under the influence of a resonant driving field. It is assumed that only one 

realization of this process is available. We show that it is possible to approximately visualize 
in real time the evolution of the system as far as it is given by |c2(i)| 2 . For this purpose we 
use a sequence of particular unsharp measurements separated in time. They are specified 
within the theory of generalized measurements in which observables are represented by 
positive operator valued measures (POVM). A realization of the unsharp measurements 
may be obtained by coupling the two-level system to a meter and performing the usual 
projection measurements on the meter only. 



& . 

§h; 1 Introduction 

Recently, considerable progress has been made in the experimental and theoretical study of 
single quantum systems in atomic traps and optical cavities as well as in the context of solid- 
state physics (coupled quantum dots). In many cases the physical situation may be described 
by means of a two level system with states |1) and |2) under the influence of a periodic driving 
potential V(t). The resulting motion of the normalized state 

|^) = Cl |l)+c 2 |2) (1) 

involves oscillations of the probabilities |ci(t)| 2 and | C2 ) | 2 called Rabi oscillations. The usual 
way to measure the dynamics of |c2(t)| 2 employs many projection measurements of an observable 
with eigenstates |1) and |2). For this purpose an ensemble is prepared in the initial state and a 
projection measurement is carried out at time to on each member of the ensemble leading to the 
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determination of | C2 (^o) 1 2 - Repeating the procedure for different times t one obtains the Rabi 
oscillations of |c2(£)| 2 . 

This procedure fails if there is only one single quantum system available and the objective is 
to visualize its motion in real time. Imagine, e.g., there is a single two-level system performing 
only once a hundred Rabi oscillations. Is there a method to visualize them? There are two 
complementary conditions such a method must meet. The measurement must not disturb the 
system too much, so that its evolution remains close to the undisturbed case. On the other 
hand, the coupling of the measurement apparatus to the two level system should nevertheless 
be so strong that the readout shows accurately enough the modified dynamics of the system 
(including the disturbance by the measurement). 

This can clearly not be achieved by standard measurements which project onto an eigenstate 
of the measured observable. They disturb the system too strongly. For instance, the well known 
Quantum Zeno effect shows that it is possible to detect perfectly the modified dynamics of a 
system at the price of loosing all information about the undisturbed motion. The measurement 
disturbs the dynamics so much as to make the system stay in one of the eigenstates of the ob- 
servable. A more promising strategy to fulfill our demands is to carry out appropriate sequences 
of generalized measurements. Generalized measurements differ from projection measurements in 
that the probabilities of their outcomes are in general given by the expectation values of positive 
operators instead of projection operators. 

The theory of generalized measurements in which observables are represented by positive 
operator valued measures (POVM) has intensively been studied For a survey see Busch et 
al. ||. These measurements are in their generic form also called unsharp, as we will do, or non- 
ideal, weak, soft or fuzzy, thereby paying attention to one of their different properties. They have 
been investigated theoretically in connection with the problem of joint non-ideal measurements 
of incompatible observables, c.f., Busch et al. 0. Recent quantum optical measurements like the 
Haroche- Ramsey set-up have been analyzed too with respect to the question of complementarity 
by de Muynck and Hendrikx || using the concept of generalized measurements. Martens and 
de Muynck Q address the question of which measurements extract the most information about 
the system. Information extraction and disturbance has been widely discussed in quantum 
information theory, as e.g., in the context of continuous measurement and feedback control by 
Doherty et. al. |J. We will deal instead with sequences of discrete unsharp measurements which 
are separated in time and show that, if appropriately chosen, they may provide a real time 
visualization of Rabi oscillations. 

In the context of continuous quantum measurements of energy the question of visualization 
of oscillations has been addressed in in a phenomenological approach, which was based on 
the application of restricted path integrals 0. In || it had been indicated, that for continuous 
measurements a correlation may exist between the time-dependent measurement readout and 
modified Rabi oscillations between energy eigenstates. This visualization of a state evolution 
has been numerically verified in ||, to our best knowledge for the first time (compare also 
Oscillations in coupled quantum dots measured by a quantum point contact are treated in an 
approach to continuous quantum measurements using stochastic master equations by Korotkov 
fLOf . For a quantum trajectory approach see Goan et al. |TT]. 



A realization scheme for the continuous fuzzy measurement of energy and the monitoring of 
a quantum transition has then been given in |l2j . The intention was to present a microphysical 
basis for the phenomenologically motivated continuous measurement scheme. Inspired by this 
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article we give below an independent treatment which is solely based on successive single-shot 
measurements of the generalized type (POVM). All concepts necessary to introduce an appro- 
priate measurement readout, to justify this choice and to specify different measurement regimes 
are based on these single measurements. This is also the case for the numerical evaluation which 
makes use of the simulation of single unsharp measurements and the otherwise undisturbed dy- 
namical evolution (Rabi oscillations) between the measurements. Our intention is to keep our 
treatment as general as possible but nevertheless readable for those who may be able to design 
experimental realizations. We will give an application to a particular quantum optical setup in 
a subsequent paper Q| . 

This paper is organized as follows: In Sect. |2| we specify the particular subclass of unsharp 
measurements on which our considerations are based. In Sect. |3| a succession of these unsharp 
measurements is studied. In addition we introduce the concept of a 'best guess' based on the 
outcome of one series of iV consecutive measurements (N-series). In Sect. § a Hamiltonian is 
introduced generating the Rabi oscillations we want to measure by means of series of unsharp 
measurements with time r between two consecutive measurements. An appropriate measurement 
readout is defined and different regimes of measurement are specified. Finally the results of the 
numerical analysis (specified in App. C) are presented in Sect. |5|. In order not to interrupt the flux 
of the presentation we postponed the detailed motivation of the particular choice of generalized 
measurements to App. A and some technical considerations concerning the measurement readout 
to App. B. 

2 Single unsharp measurements 

For later use we introduce in this section a well known special class of generalized measurements 
and discuss how they may be experimentally realized. 

2.1 A particular class of single unsharp measurements 

We discuss measurements on a two-level system (S) with orthonormed states |1) and |2), which 
are eigenstates of an observable A represented by a self-adjoint operator. This observable may 
be for example energy or spin. We perform a particular kind of generalized measurement, that 
is a measurement which can be represented by means of a Positive Operator Valued Measure 
(POVM). Its operational meaning will become clearer in Subsection |2.2| , where we describe its 
realization. The measurement is specified as follows: There are only two outcomes or pointer 
readings denoted by + and — . The duration of the measurement is St. Depending on the 
outcome the normalized initial state (^(O)) of S is tranformed due to the measurement into one 
of the related states \ip±(Sr)), 

1^(0)) \^ ± {8t)) = M ± |V(0)> (2) 

with so-called operations, which are required to be diagonal: 

M + := u+|l)(l|+it+|2>(2|, (3) 
M_ := ur|l)(l|+Mj|2)(2|. (4) 
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We add to this first requirement an important second one: uf 2 and Ui 2 are assumed to be 
positive. Here and in the following an overall phase factor exp(i</?±) is omitted without loss of 
generality. A justification for both requirements is given in App. A. Let us introduce 

ft := {utf = 1 - {u^f , p 2 : = {utf = 1 - {u 2 f . (5) 

We will characterize the measurement later on by the parameters 

Po ■= ^(Pi > Ap:=p 2 -p 1 . (6) 

According to the specifications above, we restrict ourselves to a particular class of generalized 
measurements with commuting hermitean operations, which will prove to be useful for our 
purpose. Throughout the paper we discuss selective measurements, which are non-destructive. 
We use the Schrodinger picture unless otherwise stated. Note, that the state vectors will in 
general not be normalized. 

The probabilities p± for the outcomes + and — to occur if the measurement is performed on 
the state \ip(0)) are given by the expectation value of the so-called effects 

E± = MlM ± , (7) 
E+ + E_ = 1 (8) 

according to 

P± = ( E ±)m ■ (9) 
This leads to < pi j2 < 1 and p+ + P- = 1. One may also write 

E + : = Pl |l)(l|+p 2 |2)(2|, (10) 
E_ : = (l- Pl )|l)(l| + (l-p 2 )|2)(2|. (11) 

The effects E + and £L are semi-positive operators. Because of the property (|8]) they constitute 
a POVM. It is a characteristic trait of the type of measurements in question that the effects are 
in general not projections: E± ^ E±. We call this the genuine case. 

Only in the limiting cases with p\ — \ and p 2 = (or with indices 1 and 2 interchanged), 
i.e. |Ap| = 1, we obtain a projection operator valued measure (PVM). Looking at the resulting 
states \ijj±(8r)) in this limit we see that they are identical with |1) or |2). Accordingly, if the 
initial state ^(O)) is |2), we obtain reliably + as measurement outcome and — for |1) (or with 
indices 1 and 2 interchanged). This well known von Neumann measurement is called a sharp 
measurement. 

The opposite limit is characterized by |Ap| 1. Then E± are nearly proportional to the 
identity operator and the probabilities p± of the outcomes become nearly independent of the 
initial state ^(O)). Even if the initial state is |1) (or |2)), there are finite probabilities for 
both outcomes + and — . Because in the genuine case the generalized measurement does not 
allow a fully reliable conclusion regarding the initial state, it is called an unsharp (or non-ideal) 
measurement. With E± of the form (|10"D and (|LT|) it is generally accepted to refer to it as an 
unsharp measurement of the observable A . Because it is not evident what this notation could 
mean from an operational point of view, we discuss an explicit realization below. In the limit 
|Ap| 1, the measurement changes the state only a little as can be seen from equations (|3]) 
and (f|). Regarding its influence on the state, an unsharp measurement may also be called weak 
or soft. Below we will make use of the fact that single measurements with different degrees of 
weakness may be introduced by choosing the parameters po and |Ap| appropriately. 
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2.2 Realization of unsharp measurements 



It is known that the unsharp measurements above can be realized by coupling the two-level 
systems S to a meter M and performing a projection measurement on M only (c.f. Neumark's 
theorem ||). It is to be expected, that in this way experiments with unsharp measurements of 
the class we are discussing may in practice be performed for example in the domain of quantum 
optics or condensed matter physics. We give therefore some details: 
The two-level system in the state 

= Cl |l) + c 2 |2) (12) 

interacts with another two-level system with orthonormed states |$±) which acts as the meter M. 
Its initial state is | $(0)). Hence the initial state of the compound system is |\l/(0)) = \ip(0))\§(0)). 
It is then assumed that the interaction lasting for the time 5t results in the particular unitary 
development 

|*(0)> - mSr)) (13) 

= «+ Cl |i)i$+)+wrci|i)i$_) 

+ u 2 + c 2 |2)|$ + >+u 2 c 2 |2>|$_} 

with positive uf t2 ,u^ 2 and ||^(0)|| = ||^(5r)|| = 1. 

After St a sharp measurement on the meter alone is performed which, depending on the 
outcome + or — , projects on the state |<E>+) or |$_) with P± := 1® |<&±)($±|. These two pointer 



readings + and — are taken as the measurement outcomes of Section |2TT| . The readout + occurs 
with probability 

which agrees with p + of (^j) . In this case the state of the compound system after the measurement 
is the product state 



P+ = (P+W) = «) 2 | Cl | 2 + (ut) 2 \c 2 \ 2 , (14) 



\*+(5t)) = |¥(*r)) = « Cl |l) + M2 + c 2 |2»|$ + > = |^ + (*r))|* + ) (15) 

with \ip + (5r)) as given in equations (^) with (|3|) and (^). In order to obtain the final state for 
the measurement readout — one simply has to interchange the indices + and — . The class of 
unsharp measurements of Section is thus reconstructed. 

The device specified above shows clearly that the unsharpness in question is of genuine 
quantum nature. It cannot be traced back to an imperfect measurement procedure, which itself 
may for example be caused by an imperfectly fixed pointer showing from time to time + although 
the meter state is The reason is instead that the quantum dynamics leading to the unitary 

development ( |T3"D does not correlate |1) only with |$ + ) and |2) only with |$_) but allows the 
appearance of the states |2)|$ + ) and |1)|$_) in the superposition. 



3 N-series of successive unsharp measurements 
3.1 One N-series 

We will now perform on the same single system N unsharp measurements of the type specified 
above with the same parameters in an immediate succession. This is called an N-series. Later 
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on we will allow the measurements to be separated in time by r. A motivation to use N-series 
will be given below. Since the operations M + and M_ commute, the final state after the N-series 
will not depend on the order of + and — results. For any particular sequence of results with a 
total number iV + of results + the state is transformed for r = according to 

mO))^M(N + ,NM(0)) (16) 

with 

M(N+, N) := M+ + M { _ N ~ N+) . (17) 

Now we discard the information about the order of the results, i.e., we restrict ourselves to 
the information that the total number of + results is N + , regardless when they occurred in the 
sequence. As measurement outcome attributed to the total N-series we take in the following the 
relative frequency of positive results 



In order to transcribe for this type of unsharp measurement the scheme of Section [2.1| , which 
is based on the equations (0), (|7|) and (^) for non- normalized state vectors, we have to work 
out the related operations and effects. The probability p(N + ,N) that N + positive results are 
measured in an N-series is times the probability that a particular ordered sequence of N + 
positive and N — N + negative results is obtained: 

p{N + ,N) = (E{N + ,N)) m , (19) 

where the effect E(N + , N) is given by 

E(N+, N)=(^] M t (A+, N)M{N+, N) (20) 




Pi 



^ + (l-Pi) (JV - JV+) |l)(l|+pf + (l-p 2 ) (JV - JV+) |2)(2| 



The corresponding probability is 

P(N + ,N) = ^(P^il-P.r^hl 2 (21) 

+ ^(l- P2 )(^+)|c 2 | 2 ) . 
To complete the scheme we write the pure state resulting after the N-series in the form 



' N 



\M(N + ,N)\i>(0))- (22) 



In order to quantify by how much an N-series of measurements affects the state, we calculate 
for arbitrary values of p\ and p 2 the fidelity, which shows how closely the post-measurement 
state resembles the pre-measurement state. The fidelity F(N) between the pure state ^(O)) 
before the measurements and the state 

N 

P= E l^ + )(^AT + | (23) 
N+=0 
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after the N-series is equal to the square of root of the overlap between \ip(0)) and p 



F(J\0 == VMHVO ■ ( 24 ) 
Please note, that \ipN+) is unnormalized. After some algebra one finds 



F(N) = ^l-2\ Cl \ 2 \c 2 \ 2 (l-b), (25) 

where 



N 



y/pm + y/(l-pi)(l-P2)) ■ (26) 

Equation fl25l) shows that for all choices of the parameters pi and p 2 the fidelity assumes its 
minimum for |ci| 2 = | C2 1 2 = 1/2 and its maximum for |ci| 2 = 1 and | C2 1 2 = (or with indices 
interchanged). 

Fidelity may serve as a direct measure for the weakness of the influence of an N-series of weak 
measurements on the pre-measurement state. This becomes evident by looking at the limiting 
cases: for extremely sharp measurements with pi = 1 and p 2 = (or vice versa) the fidelity 



equals the fidelity F = — 2\c 1 \ 2 \c 2 \ 2 of a projection measurement. The maximal fidelity 
F = 1 is obtained for infinitely weak measurements with pi = p 2 . 



For an infinite N-series the limit limjv^oo F = yl — 2|c 1 | 2 |c 2 | 2 equals for all values of pi = 1 
and P2 = with p% ^ p 2 the fidelity of an projection measurement. That such an iteration of 
measurements produces an eigenstate has theoretically been discussed in [|15|| . For an experimen- 
tal scheme which realizes a projection measurement by an iteration of unsharp measurements 

m. 



sec 



3.2 Best guess for the outcome of one N-series 

For the moment we refer again to repeated measurements on the same initial state \ip(0)) (en- 
semble approach). The statistical expectation value of r = N + /N reads 

N 

£(r) = £ r(N + )p(N + ,N) (27) 

Af+=0 

= Pl\ci\ 2 + | C 2 1 2 • 

The latter equation follows with (pl|). In a similar way we obtain the variance a 2 (r): 

o 2 (r) = | Cl | 2 |c 2 | 2 (Ap) 2 + f| Cl |V(l -Pi) + |c 2 | 2 p 2 (l -P2)) • (28) 



N 

Based on (p7[) it is easy to relate \c 2 \ 2 of the initial state ^(O)) with the expectation value £{r) 
of an N-series of unsharp measurements starting with |-0(O)}: 



o £(r)-pi 



Please note, that both quantities \c 2 \ 2 and £ (r) refer to an ensemble represented by the state 
|V(0)>. 
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If only one N-series is measured as will be the case below, only a best guess G 2 for the quantity 
|c 2 | 2 can be obtained. Equation fl2"E| ) suggests to choose it as 

G 2 :=^. (30) 

The standard deviation of the best guess G 2 is then given by 



<T(G 2 ) = ^Pr (31) 



IAp| 



|ci| 2 pi(l - pi) + |c 2 | 2 p 2 (l - P2) 



Y 1 11 n (A P y 

It specifies how accurate |c 2 | 2 of the pre-measurement state ^(O)) can be estimated by means 
of one N-series. cr(G 2 ) can be decreased by increasing iV or |Ap|. 

In case of a single measurement (N = 1), the standard deviation <r(G 2 ) assumes its minimal 
value |ci||c 2 | for |Ap| = 1, i.e., for a projection measurement (pi = 1 and p 2 = or vice versa). 
For all other choices of the parameters p± and p 2 (|Ap| < 1) and other choices of N, there is 
a larger standard deviation. This justifies to call these measurements unsharp, because they 
provide less reliable information about |c 2 | 2 . 

As we will see below, there are physical situations in which the parameters are fixed and have 
to be chosen with |Ap| < 1. Eq. (|3l| ) shows that in this case the standard deviation decreases 



with increasing number of repetitions N so that it is favorable to refer to N-series. This is the 
reason why we will now make use of them. 

4 Measurement of a dynamically driven state by a se- 
quence of N-series 

4.1 Best guess as measurement readout 

We refer in the following not to an ensemble but to one single two-level quantum system. It is 
driven by an external periodical influence with frequency uj to perform resonant Rabi oscillations, 
if it is not disturbed by measurements. In the Schrodinger picture its undisturbed dynamics is 
therefore given by 

H = ai|l)(l| + a 2 |2)(2| + — — ( 1 2> < X | exp{-iujt} + |1)(2| exp{iujt}) (32) 

with Q R = 2n/T R and Rabi period T R . Because of the resonance we have: to = (a 2 — ai)/h. 

It is our intention to visualize the undisturbed behaviour in time of |c 2 (t)| 2 , i.e., in our case 
the Rabi oscillations. A sharp measurement implies a drastic change of IV^o)) into |1) or |2). 
A succession of projection measurements will in fact result approximately in the Zeno effect. 
Unsharp measurements can be weak and may change the state less. If their parameters are 
properly adjusted, one can therefore expect that these measurements are superior for our pur- 



pose. Accordingly many unsharp measurements of the type specified in Sect. 2.1 are performed 
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on this system with time r between two consecutive measurements. The unitary development 
between the measurements is determined by the Hamiltonian H . It is assumed that the time 
t is much larger than the duration 5t of a single unsharp measurement (r ^> St). Again we 
bundle up N unsharp measurements to an N-series of total duration At := Nt, so that in fact a 
sequence of N-series is obtained with outcomes r(t m ) at t m = mAt with m = 1 , 2 . . .. In order 
to resolve the Rabi oscillations we require At <C T R . As outlined in Appendix B the best guess 
G 2 {t m ) related to r(t m ) by 

G 2 (t m ) = r -MzP± (33) 

can still be used as approximation for |c2(t m )| 2 even if a driving influence is present. However, 
it will only be a good approximation if the number N of measurements within a N-series is not 
too large. More precisely, N has to fulfill inequality ([45]). The sequence {G 2 {t m )) obtained from 
the outcomes of the successive N-series serves as the measurement readout. 



4.2 Level resolution time and the different regimes of measurement 

How unsharp (or weak) must the underlying single measurement of Section ^| for our purpose 
be? A measure for the strength of the driving influence in fl32|) is given by the frequency Qr. A 
strong influence corresponds to a short Rabi period of the oscillations between the two levels |1) 
and 1 2). On the other hand measurements on the system, even if they are unsharp, modify and 
hinder this primary time development of the state of the system. We have therefore two contrary 
influences on the two-level system. For a comparison we want to introduce a characteristic time 
which represents the reciprocal strength of the disturbing influence of the measurement alone. 
To this end we discuss in a first step again one N-series with r = 0. 

The two states |1) and |2) of the two-level system correspond to | C2 1 2 = and | C2 1 2 = 1- 
According to ( |3"1"D the standard deviation a(G 2 ) of the best guess G 2 improves with increasing 
N. How large must N be in order to obtain a standard deviation smaller than 1, so that one 
can discriminate between the two values of \c 2 \ 2 above? We have to demand 

1 > o-\G 2 ) . (34) 

a{G 2 ) depends on \c\\ and |c 2 | of the initial state. But this requirement should be fulfilled for 
any initial state. We may therefore take in (|34]) the maximum value of <j(G 2 ) 

"4 + iV 2(KpJ ' [6b) 

which amounts to 

In fact we are interested in unsharp measurements within the N-series which are repeated after 
a time r > 0. In order to discuss the influence of the measurements separately, we consider the 
absence of the driving influence (g — 0). Then the operator of the dynamical evolution between 
the measurements commute with M±. Therefore independent from the initial state altogether a 
finite interval of time not larger than 

T ir = rN lr = r 3{Ap)2 (37) 
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is needed in order to obtain a variance of the related best guess so small as to resolve the levels 
with high probability. We call T\ r the level resolution time. 

We return to the situation when in addition the driving external influence is active too. The 
complete setup is then characterized by two time scales: the level resolution time T; r and the 
Rabi period Tr. There are three physically different regimes: for T\ r <C Tr the succession of 
measurements is strong or sharp enough to push the state to one of the states |1) or |2) before the 
Rabi oscillation becomes manifest. Connected with this is a small variance. Therefore G2(t m ) 
will correspond well to the disturbed |c2(t m )| 2 but the readout will not reflect the underlying Rabi 
oscillations. We call this the quantum jump regime for reasons to be seen later. For Tr <C T\ t 
on the other hand, the measurement influence on the state dynamics is so weak, that the state 
and therefore |c2(t m )| 2 shows the Rabi oscillations. But because of the large variance of G2, 
a single readout G^m) will in general not be significantly related to |c2(t m )| 2 but will show a 
fuzzy behaviour. The state motion thus cannot be read off from G^t^m)- We call this the Rabi 
regime. Of particular interest is the intermediate regime characterized by Ti t /Tr m 1. One can 
expect that this regime shows two features: on one hand the state motion is still close to the 
Rabi oscillations and on the other hand the readout G^ifm) follows closely the |c2(t m )| 2 of the 
state. In this case the state motion would be made visible as intended, without modifying too 
much the original Rabi oscillations by the influence of the measurement. 

The numerical simulations show that the intermediate regime can rather be characterized by 
the condition 

/ := 3nT lr /T R « 1 (38) 
than by T^/Tr « 1. In the following we use / instead of TJ r as measure of fuzziness. 

5 Results of the numerical analysis 

We discuss the three regimes by considering one typical example for each regime. The charac- 
teristic traits can be read off from Figs. ^ [3[ In all three cases we have chosen r = 0.002Tr and 
N = 25, so that condition (^5| ) is fulfilled. Figs, pi- ||a show as function of t/T R 

• the |c2| 2 -curve of the state motion (black), 

• the GVcurve representing the measurement readout (gray) and 

• the processed GVcurve obtained after noise reduction (dashed). 

The noise reduction is sketched in Appendix C. The related Figs. 0b- |3|b show the respective 
Fourier spectra gained by discrete Fourier transform of the | C2 1 2 -curve (black) and the measure- 
ment readout G 2 (grey). The peaks at k = are caused by the non- vanishing average value of 
the curves. In the three regimes we find characteristically different results: 

5.1 Quantum jump regime 

With fuzziness / = 0.07 there is a rather strong influence of the measurements on the time 
evolution of the state. The single measurement is still rather sharp and close to a projection 
measurement. Accordingly, we observe in Fig. [TJ, that the Rabi oscillations of | C2 1 2 , which would 
appear in the absence of the measurements are completely destroyed. The |c2| 2 -curve (black) 
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shows jumps between extended periods with values approximately equal to or 1. This means 
that if the state of the system is close to an eigenstate |1) or |2), it has a high probability to 
stay there. The GVcurve of the measurement readout (gray) reflects very well the |c2| 2 -curve of 
the state. The peaks of the |c2| 2 -curve pointing upwards and downwards follow the direction of 
the corresponding peaks of the G 2 -curve. The processed G 2 -curve (dashed) reflects the quantum 
jump regime even better. In this regime there is a rather good agreement between the readout 
and the |c2| 2 -curve, but there is no characteristic indication of the undisturbed dynamics. The 
Fourier spectra show a very good correlation, as can be seen in Fig. |I|b. 

5.2 Rabi regime 

The other extreme can be found in the Rabi regime. We have chosen a high fuzziness / = 63. The 
single measurements are so unsharp and their influence on the motion of the state so weak, that 
the Rabi oscillations remain nearly undisturbed. This can be seen from the |c2| 2 -curve (black) in 
Fig. ^|a. In contrast to this the GVcurve of the measurement readout (grey) is very fuzzy. The 
single measurement results lie far away from the interval [0, 1]. The processed GVcurve (dashed) 
shows also no correlation, neither in phase nor amplitude with the |c2| 2 -curve. Turning to the 
Fourier spectra (Fig. ^]b), we see that the Rabi frequency of the |c2| 2 -curve (black) is represented 
well. On the other hand, no frequency can be attributed to the measurement readout (gray 
curve). In this regime the measurement readout does not allow a conclusion on the motion of 
the state as represented by the |c 2 | 2 -curve. 

5.3 Intermediate regime 

Finally we turn with a fuzziness of / = 0.98 and |Ap| = 0.08 to a parameter choice, which 
leads to representative curves for the intermediate regime. According to the Fourier spectrum 
(black) in Fig. |3|b the |c2| 2 -curve is dominated by one frequency, which differs only little from 
the Rabi frequency Tr. In this regime the |c2| 2 -curve, which would characterize the state motion 
if no measurement were performed, is not much modified. The important point is that now the 
GVcurve of the measurement readout (grey) shows a correlation with the |c2| 2 -curve regarding 
the phase: if the |c2| 2 -curve has a maximum, the G^-curve is in the upper domain and for a 
minimum correspondingly. Regarding the Fourier spectrum of the measurement readout (grey 
curve in Fig. ^|b), there is still a lot of noise, but a clear peak indicates the frequency of the | C2 1 2 - 
curve. The measurement readout follows therefore closely the |c2| 2 -curve of the state motion in 
phase and frequency. This is even more convincing for the processed G^-curve (dashed curve in 
Fig. |a). 

6 Discussion and conclusions 

The normalized state \ip{t)) = Ci(£)|l) +C2(t)|2) of a single two-level system performs oscillations 
under the influence of a resonant time dependent driving field. It is assumed that there is only 
one realization of this process available. We ask the question whether it is possible to visualize 
in real time the evolution of the system as far as it is given by |c2(t)| 2 . 

This result can clearly not be obtained by means of standard measurements which result in 
projections. Generalized measurements, can be discussed by means of positive operator valued 
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measures (POVM). For our goal we restrict to a particular subclass of unsharp measurements, 
which have a simple structure so that an experimental realization should be possible. The pointer 
values of this measurements are + and — . We perform a sequence of measurements separated in 

time leading to the pointer readings H 1 !- + ■■■■ It can be shown that it is favorable to 

transform a series of N pointer readings to one point of a measurement readout which is defined 
as best guess G 2 (t) for |c 2 (t)| 2 . 

The single unsharp measurement must be weak enough as not to disturb the evolution of 
the system too much. On the other hand it must not be too weak in order to give information 
about the state of the system. We determined analytically the domains for the different open 
parameters such that this 'intermediate' regime of measurement can be reached. The numerical 
analysis finally shows that for this intermediate regime the real time visualization of |c2(i)| 2 
can be obtained to a satisfactory amount. The readout (^(i) reflects rather well the resulting 
I C2 ) | 2 , which itself shows only a slight modification of the undisturbed Rabi oscillations. 

Note that above the Rabi frequency of the unitary development was assumed to be known. 
If this is not the case a feed-back mechanism must be introduced so that after a certain time 
from the beginning of the measurements a visualization can be reached. This project is subject 
of future research. 
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Appendix A: How to choose the operations 

In order to detect the dynamics of | C2 1 2 we choose measurements with particular operations 
corresponding to a particular POVM. We justify this choice as follows: 

Since our measurement scheme contains consecutive measurements, it is important to con- 
sider in which state the system is left after one measurement. We require the operations, which 
describe the transformation of the state up to a phase factor, to be such that the states |1) and 
\2) are not changed due to the measurement. Hence these states should be eigenstates of the 
operations M + and M_ in agreement with (|3|) and (Q). 

We now address the question why M± should be positive operators, thus it^2 ■> ^i 2 ^* i n 
(Q) and (second requirement of Sect. |2~T| ). In the context of quantum feedback control one 
makes use of the fact that the operations may be decomposed into 'modulus' and 'phase', in our 
case: 

M± = U±\M±\, (39) 



where U± are unitary operators and \M±\ := y M±M±. While the 'modulus' |M±|, because 
of its relation to the effects E± = \M±\ 2 , is connected to the acquisition of information, the 
unitary operators U± represent rather an additional Hamiltonian evolution, not leading to any 
increase of information about the system (compare Wiseman |T7| and Doherty et al. pfl). If 



thus U± cannot be used as a feedback to compensate the influence of M± , it only adds to the 
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disturbance by the measurement and should be set U± — 1 (a global phase factor exp(i^ ± ) is 
always neglected). 

That the additional evolution represented by U± is indeed obstructive in our case can be 
made clear as follows. To be diagonalizable the M± have to be normal, i.e., [M±, M±] = which 
amounts to [|M ± | 2 , U±] = 0. Hence U± and |M±| have a simultaneous basis of eigenvectors, 
which has to be |1) and |2), because M± is required to be diagonal in that basis. Now if U± is 
diagonal with respect to |1) and |2), it can be written as: 

U ± = exp(~a z 9 ± ) (40) 

for some angle 9 and a z as Pauli operator. In order to illustrate the corresponding different 
influences on the state motion, we refer to the Bloch sphere. The resonant Rabi oscillations 
which our system performs in the absence of measurements can be represented by a rotation 
of the normalized Bloch vector in the y-z-plane (given that the initial state lies in this plane) 
. The measurement induces a change of the state given by M±. The contribution due to \M±\ 
leads to a jump of the Bloch vector towards one of the 'poles' on the Bloch sphere at (0, 0, 1) 
and at (0,0,-1). The Bloch vector thereby does not leave the plane spanned by the state 
before the measurement and the z-axis. U± on the other hand forces the Bloch vector to leave 
this plane. As a consequence U± cannot compensate the change due to \M±\ and the influence 
of \M±\ is inevitable. Beyond that, a unitary evolution U± ^ 1 disturbs the original Rabi 
oscillations heavily in a succession of measurements by leading the Bloch vector out of the y-z- 
plane. Therefore we have to choose U± = 1 in order to disturb the Rabi oscillations by the action 
of the operations M± as little as possible. In an experimental situation with nontrivial unitary 
transformations U± the state change due to U± could be compensated by means of feedback 



(compare Wiseman [17 



Appendix B: Upper limit for N 

In this Appendix we describe under which circumstances the concept of the best guess as intro- 
duced in Sect. [5] can also be used for N-series with dynamical evolution due to the Hamiltonian 
in (|32"D between the single measurements. 

Consider a N-series of unsharp measurements which are separated by the unitary evolution 
according to (|32|) . In the interaction picture the state after the first N-series thus reads: 

n (Mfcexp^-^r))^)) (41) 

with Mfc = M + , M_ and Hj := M7r(|2)(1| + |l)(2|)/2. We have to discuss for this new situation 
again the concept of a best guess. If we may commutate unitary evolution and operations 
representing the change due to the measurement, so that 

\ip(At)) « exp (—HjAtj M(N +} N)\ip(0)) (42) 



with 

' N 



M(N+,N) 



\M(N + ,N), (43) 
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we are back to the situation described in Section |3]^ and the best guess G 2 (t m ) of (^) may be 
taken as a good candidate to approximate \c2(t m )\ 2 . The underlying commutators are 



M. 



± , exp 



a y {uf 



,7TT, 



sin 



A' 



(44) 



where represents the Pauli operator. Although the commutators are small for t/Tr <C 1, 
their contribution accumulates and is no longer negligible if the number N of measurements 
within the N-series is too large. 

A lengthy calculation shows that for Nt/Tr <C 1 the term (IT) may be approximately written 
in the form (f42|), if 



(N-iy < 



(45) 



2max{|-U2 ~~ u i |> \ u 2 ~ u i 1} 7lT 

Here it was assumed that the matrix elements of the operation M + fulfill Ui > u\ ■ Otherwise 
the indices 1 and 2 in ( f4"5| ) have to be permutated. 



Appendix C: Simulation procedure and noise reduction 

In this appendix the procedure is described we employed to simulate the measurements of Rabi 
oscillations on the computer. It is denoted here in form of four instructions: 

1) Evolve initial state unitarily according to the undisturbed dynamics of the system; \i/j(0)) — > 

|V(t)) = exp{-!if 7 T}|^(0)) . 

2) Simulate single unsharp measurement in two steps. First generate a random outcome by 
a Monte Carlo method according to the probabilities p± = {E±)^ T y Then change the state 
depending on the outcome — > M±\ip{r)). 

3) N-series: Carry out instructions 1) and 2) N times, every time replacing the initial state by the 
state obtained after the previous unsharp measurement. After the N-series record | C2 (ti = Nt)\ 2 
of the current state and the measurement readout, i.e., the best guess G^^i) for | C2 (^i ) | 2 - The 
best guess is calculated from the number of +- results (c.f. equation fl3*0|)). 

4) Simulate M times an N-series of unsharp measurements, always inserting as initial state 
the state obtained from the previous N-series. The total record of the simulation will then 
consist of a sequence of readouts G 2 (£ m ) and squared components of the state |c(t m )| 2 with 
t m = mAt , m = 1...M. 

The readouts from the simulated measurements can now be processed to eliminate noise. We 
used the following method. 

The sequence (G2(t m )) is expressed by means of discrete Fourier transform as G2(t m ) = J2f^o 1 a (^i) exp(ia;;t ri 
with ui = 2nl/T. From the power-spectrum (|a(a;/)| 2 ) the main peak (if present) is identified 
with the modified Rabi frequency and the noise is estimated. 

The noise is reduced by transforming the Fourier coefficients employing a Wiener filter <pi 
TB|: a{uji) — > a{uj{)(j)i. Further noise reduction is achieved by omitting higher frequencies in the 



Fourier series which cannot be related to the original Rabi oscillation. We truncated the Fourier 
series after the 2Z-th term, where 07 is the frequency of the main peak of the power spectrum. 

The resulting filtered and truncated Fourier series leads to a new sequence of best guesses 
(G^tm)) in the text referred to as 'processed measurement readout' or 'processed G^-curve'. 
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Figure 1: Quantum jump regime with fuzziness / = 0.07 and parameters Ap = — 0.3, p = 
0.5 , t = 0.002Tr , iV = 25. In Figure (a) black represents the |c 2 | 2 -curve, gray the measurement 
readout Gzif) and the dashed curve is obtained from the GVcurve by noise reduction. On the 
time axis multiples of the Rabi period Tr are indicated. Figure (b) shows the power spectra of 
the |c2| 2 -curve (black) and the measurement readout (gray). The frequency u is given in units 
of the Rabi period Qr. Both spectra agree very well. 
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|a(co)| 




Figure 2: Rabi regime with fuzziness / = 62.8 and Ap = 0.01 ,p = 0.5 , r = 0.002T R , N = 25. 
In Figure (a) the numerically obtained |c2| 2 -curve (black) agrees with the undisturbed Rabi 
oscillations. On the other hand | C2 (t) | 2 is neither significantly correlated with the measurement 
readout (G^-curve, plotted in gray) nor with the processed GVcurve (dashed) obtained from the 
readout by noise reduction. Fig. (b) displays the power spectrum of | C2 ) | 2 (black) and the power 
spectrum of the measurement readout (gray). For graphical reasons the power spectrum of the 
measurement readout has been multiplied by a factor of 1/5. Both spectra are not correlated. 
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Figure 3: Intermediate regime with fuzziness / = 0.98 and Ap = 0.08, p = 0.5, r = 0.002Tr, 
N = 25. Fig. (a): There is a high correlation between the |c2| 2 -curve (black) and the measure- 
ment readout G2(t) (gray) leading to a very good correlation after the noice of the readout has 
been reduced (dashed curve). The |c2| 2 -curve indicates that the motion of the state has been 
disturbed by the measurement. But the Rabi oscillations are still recognizable. Fig. (b): The 
main peak of the power spectrum of the |c 2 | 2 -curve (black) coincides with the main peak of 
the power spectrum of the GVcurve and indicates approximately the Rabi frequency Qr of the 
undisturbed motion. 
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